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ABSTRACT 

This  report  documents  work  initiated  during  the  first  nine  months  of  a 
jor  new  five-year  program  which  is  focused  cu  some  of  the  fundamental  problems 
in  ferroelectric  and  ferroelastic  crystals  which  are  of  special  relevance  to 
the  design  and  development  of  new  and  improved  materials  for  electromechanical 
transducers . 

Electrostriction  is  the  basic  coupling  mechanism  between  dielectric  and 
mechanical  properties  in  all  present  ceramic  piezoelectric  transducers,  yet  there 
is  no  adequate  theoretical  description  of  this  phenomenon,  and  reliable  experi- 
mental values  are  not  available  for  the  electrostriction  constants  of  many 
simple  insulators. 

For  experimental  measurement  of  electrostr ictive  deformations  an  ultra- 
sensitive AC  dilatometer  is  being  developed.  The  present  instrument  will  resolve 
Hi  -mui  ^[0 .005a|  at  7.5  Hz,  and  is  presently  being  used 


displacements  of  5^0  ^nwtere 

to  measure  electrostriction  constants  for  a range  of  alkali  halide  crystals. 

A comparative  evaluation  of  lattice  theoretical  models  previously  used 
to  explain  the  electrostriction  constants  in  alkali  halides  shows  the  inadequacy 
of  all  current  models  and  of  the  basic  lattice  dynamical  framework  used.  In 
order  to  develop  a more  adequate  approach  the  lattice  theory  of  the  elastic 
dielectric  was  begun  for  an  anharmonic  shell  model,  in  which  the  nonlinearity 
of  the  electronic  polarizability  is  included  in  a consistent  manner  in  addition 
to  the  anharmonicity  of  the  short-range  overlap  and  of  the  coulomb  forces.  If 
successful  for  alkali  halides,  the  theory  will  be  applied  to  piezoelectric  and 
ferroelectric  prototypes  of  simple  crystal  structure. 

For  the  important  practical  electrostrictor  lead  magnesium  niobate  (PMN) 


single  crystals  have  been  grown  which  are  large  enough  to  permit  measurements 
of  the  electrostriction  constants.  A new  digital  program-controlled  flux  growth 


system  has  been  Installed  and  is  being  used  to  grow  single  crystals  of  composi- 
tions in  the  PbZrO^PbTiO..  solid  solution  family. 

New  replication  techniques  are  being  developed  for  the  synthesis  of  ceramic: 
polymer  composites  with  controlled  interphase  connectivity  designed  to  optimize 
the  properties  for  piezoelectric  voltage  mode  applications.  Synthesis  tech- 
niques have  also  been  developed  to  generate  ferroelectric  powders  in  the  bismuth 
oxide  layer  structure  compounds  which  have  very  marked  shape  anisotropy,  so  as 
to  permit  the  fabrication  of  grain  oriented  piezoceramics  without  the  expensive 
and  time  consuming  step  of  uniaxial  hot  pressing. 

Work  has  also  been  initiated  on  the  growth  of  bicrystals  of  ferroelectric 
triglycine  sulfate  and  of  paraelectric  potassium  tantalate  niobate  to  explore 
the  interface  phenomena  at  an  artificial  single  grain  boundary  of  controlled 
orientation.  Both  intergrown  and  diffusion  bonded  interfaces  are  under  study. 

ADAGE  computer  graphics  techniques  are  being  applied  to  the  development  of 
more  comprehensive  phenomenological  descriptions  of  f erroelectricity  in  perovskite 
tungsten,  bronze  and  bismuth  oxide  layer  structure  compounds  and  solid  solutions. 


1.  INTRODUCTION 

This  report  describes  work  carried  out  over  the  period  from  April  I to 
December  31  on  ONR  contract  N00014-78-C-0291  which  initiates  a new  five-year 
study  of  several  of  the  more  important  basic  and  applied  problems  associated 
with  the  development  of  new  and  improved  materials  and  new  understanding  of 
the  basic  phenomena  which  are  important  in  piezoelectric  transducer  systems. 

Of  necessity,  much  of  the  time  in  this  first  nine  months  lias  been  taken 
up  assembling  the  people  and  equipment  necessary  for  the  new  research  topics 
which  are  being  undertaken.  In  spite  of  this  necessary  focus,  however,  several 
important  advances  have  been  made. 

For  the  experimental  determination  of  electrostriction  in  simple  low  per- 
mittivity solids  a new  type  of  AC  capacitance  dilatometer  was  proposed.  Over 
the  past  six  months,  this  instrument  has  been  designed,  constructed,  and  is 
now  being  tested.  The  first  results  already  demonstrate  the  soundness  of  the 
original  design  principles.  The  present  sensitivity  permits  the  resolution 

o -13 

of  a displacement  amplitude  of  0.003A  (5.10  meters)  at  a frequency  of  7.5  Hz, 
thus  it  is  already  very  close  to  the  design  target,  and  adequate  for  the  pre- 
cise electrostriction  measurements  in  many  of  the  insulator  crystals  of  most 
theoretical  interest.  A paper  describing  this  instrument  will  be  presented 
at  the  32nd  Annual  Frequency  Control  Symposium  in  Atlantic  City  in  June,  1979. 

The  theoretical  studies  of  electrostriction  are  also  well  underway.  Dr. 
Achar,  working  in  close  cooperation  with  Professor  Barsch,  has  recalculated 
the  electrostriction  constants  for  alkali  halides,  using  the  three  presently 
available  theoretical  models  but  with  the  most  recent  values  of  the  input  para- 
meters. These  calculations  clearly  exhibit  the  poor  agreement  from  presently 
available  theories,  and  confirm  the  disturbing  feature  that  the  more  sophisticated 
models  do  not  represent  significant  improvement  over  the  earliest  point  charge 
treatments . 
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A feature  of  the  earlier  model  treatments  which  is  believed  to  be  most 
important  is  the  neglect  of  the  shell  polarizability  in  the  shell  models  applied. 
The  theory  is  now  being  reformulated  to  take  this  into  account.  Phenomenologically 
it  is  evident  that  in  the  systems  which  have  been  measured,  low  permittivity 


solids  generally  have  higher  Q constants  and  high  permittivity  solids  have 

lower  values,  where  the  Q..,.  are  the  polarization  related  electrostriction  cons- 

ljkl 

tants.  Since  there  is  more  than  2 orders  of  magnitude  in  the  range  of  Q values, 


it  is  obviously  very  important  to  understand  whether  this  coupling  is  fundamental. 


If  higher  K materials  could  be  found  which  have  very  high  Q...  , values,  their 

i]kl 

application  would  make  a major  improvement  in  transducers. 

Work  on  two-phase  ceramic : plast ic  and  ceramic : ceramic  composites  is  being 
carried  forward  on  two  parallel  approaches.  Dr.  Brendan  Hardiman,  who  lias  joined 
the  group  from  A.G.  Telefunken  Laboratories,  is  exploring  new  techniques  for 
replicating  3:1  connected  ceramic : plast ic  composites  in  which  the  scale  of  the 
macrostructure  can  be  tailored  over  a very  wide  range.  The  method  uses  a lost- 
wax  casting  process,  but  now  with  a reusable  template  which  can  be  controlled 
very  precisely.  In  an  alternative  strategy,  we  are  exploring  extrusion  tech- 
niques for  making  ultrathin  PZT  rods  (100  to  200  pmeters)  which  can  then  be 
fabricated  in  the  flexible  green  state  into  many  of  the  complex  interconnected 
structures  required. 

The  study  of  grain  oriented  ceramics  has  focused  on  the  task  of  generating 

fine  powders  of  the  bismuth  oxide  layer  structure  compositions  whicli  have  very 

marked  shape  anisotropy.  By  reacting  the  mixed  oxides  in  a eutectic  KCl:NaCl 

salt  bath,  we  have  found  that  Bi.Ti.O,  , and  Bi.,WO,  can  be  produced  as  a micron- 

4312  2b 

sized  powder  with  very  marked  pjatey  character.  Thus  without  the  expense  ot 
the  uniaxial  fiot  pressing:hot  extrusion  techniques  which  are  being  used  in 
Japan  we  can  produce  well  oriented  grain  texture. 
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To  permit  better  control  in  the  crystal  growth  studies,  a new  digital  M 
data  track  temperature  programming  system  has  been  installed  on  the  flux  growth 
furnace.  Crystals  of  Pb^MgNb.-,0^  suitable  for  the  measurement  of  the  single 
crystal  electrostrict ion  parameters  have  been  grown,  and  work  is  now  progressing 
on  the  growth  of  crystals  in  the  PZT  solid  solution  family. 

Studies  of  the  "model  grain  boundary"  have  been  initiated  using  water- 
soluble  triglycine  sulfate  and  paraeleotric  KTN  as  the  material  systems.  Bi- 
crystals of  TGS  have  been  produced  which  show  limited  areas  of  good  bonding, 
and  the  directions  for  improving  the  technique  are  being  explored.  A miniature 
hot-pressing  rig  has  been  built  for  controlled  diffusion  bonding  of  inorganic 
ferroelectrics.  Preliminary  tests  on  quartz  crystals  have  shown  the  possibility 
of  producing  good  high-strength  bonding. 

The  thermodynamic  phenomenology  for  pure  lead  zirconate : lead  titanate  solid 
solutions  developed  under  ONR  contract  N00014-76-C-051.5  was  taken  on  that  program 
to  the  stage  where  the  full  set  of  constants  up  to  sixth  rank  in  the  Devonshire 
formalism  had  been  evaluated  for  all  single  cell  compositions,  and  gave  a 
coherent  and  unified  description  of  the  dielectric,  elastic,  and  thermal  proper- 
ties . 

On  this  program,  similar  methods  using  the  ADAGE  graphics  system  are  being 
extended  to  the  practical  doped  PZTs,  and  the  PLZT  family.  The  programs  are 
also  being  reformulated  to  the  4/mmm  prototype  symmetry  so  as  to  be  able  to 
describe  the  very  large  families  of  tungsten  bronze  and  bismuth  oxide  layer 
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structure  ferroelectrics. 
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2.  ELECTROSTR I CTION 
2.1  J ntroduction 

There  is  probably  little  need  to  reiterate  at  this  stage  the  very  basic 
importance  of  electrostr ict ion  to  electromechanical  systems.  For  all  ferro- 
electrics  derived  from  centric  prototypes,  a group  which  includes  all  present 
generation  piezoceramics,  the  single  domain  piezoelectric  constants  can  be 
directly  traced  to  electrostriction  of  the  prototype.  Even  the  domain  and 
phase  boundary  contributions  to  strain  in  the  more  complex  ceramics  can  be 
traced  back  to  what  are  basically  electrostrict ive  deformations,  yet  we  have 
no  viable  theoretical  description  of  electrostriction  even  in  simple  solids  and 
almost  no  reliable  experimental  values  on  which  to  base  such  theories. 

To  rectify  this  deplorable  situation,  a two-pronged  attack  is  clearly 
necessary.  First,  the  background  of  reliable  experimental  values  for  electro- 
striction constants  of  simple  solids  must  be  built  up,  and  second  the  theoretical 
framework  must  be  established  to  give  a properly  based  atomistic  description  of 
the  phenomenon. 

In  the  experimental  area  we  proposed  to  set  up  to  measure  electrostriction 
both  by  the  direct  method  (measuring  the  change  of  sample  dimensions  under  E 
field)  and  by  the  converse  method  (studying  the  change  of  dielectric  permittivity 
under  mechanical  stress).  Over  the  past  six  months,  emphasis  has  been  on  devel- 
oping the  exceedingly  sensitive  AC  dilatometer  required  for  the  direct  study  of 
electrostr ictive  deformation. 

The  theoretical  approach  has  been  first  to  reexamine  the  earlier  theories 
of  Grindlay  and  Wong  (1969),  of  Srinivasan  and  Srinivasan  (1972)  and  of  Goyal 
et  al.  (1977)  for  electrostriction  in  the  alkali  halides.  To  recalculate  for 


representative  crystals  on  each  theory  using  a common  most  recent  set  ol  input 
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parameters  so  as  to  be  able  to  compare  properly  these  earlier  studies.  Unfor- 
tunately, this  work  has  shown  that  agreement  between  theory  and  experiment 
is  not  radically  improved,  nor  are  the  more  recent  sophisticated  theories  any 
significant  improvement  over  the  earliest  simple  point  charge  model. 

2 . 2 Development  of  the  Experimental  Technique  for  Direc t Measurement  of 
Electrostriction  in  Simple  Solids 
2.2.1  Design  of  the  AC  Dilatometer 

t 

The  dilatometer  proposed  war  based  on  a design  by  Bohaty  and  HaussuhK.1977,) 
for  an  AC  capacitor  type  dilatometer  (Fig.  2.1).  In  their  system  the  electro- 
strictive  crystal  of  interest  (a)  cut  and  oriented  in  the  correct  manner  and 
plated  with  evaporated  metal  electrodes  is  mechanically  connected  through  an 
insulating  perspex  support  to  the  upper  plate  of  the  strain  sensing  capacitor 
(b) . A standard  piezoelectric  crystal  (d)  is  electroded  and  mounted  in  a manner 
so  that  the  piezoelectr ically  generated  strain  is  communicated  in  a similar 
manner  to  the  lower  plate  of  the  strain  sensing  capacitor  (c) . If  crystal  (a) 
is  driven  with  a large  AC  voltage  at  a frequency  of  80  Hz,  plate  (b)  will  be 
excited  mechanically  with  a component  of  its  motion  at  160  Hz  (2  uj)  through 

2 

the  electrostrictive  strain  generated  in  (a)  which  will  be  proportional  to  V . 

Crystal  (d) , the  piezoelectric  standard,  is  driven  through  a frequency  doubling 
circuit  phase  locked  to  the  generator  supplying  (a) . The  piezoelectric  strain 
in  (d)  will  also  modulate  the  position  of  the  lower  plate  (c)  of  the  sensing 

1 

capacitor  at  a frequency  of  160  Hz,  in  phase  with  the  vibration  of  plate  (b) 
generated  by  crystal  (a).  When  the  phase  and  amplitudes  of  the  fields  at  (a) 
and  (d)  are  suitably  chosen,  the  two  plates  (b)  and  (c)  will  vibrate  in  exact 
synchronism  producing  a net  zero  modulation  of  the  separation  and  thus  of  the 


capacitance  of  be  at  160  Hz. 


BRASS 


AC  differential  dilatometer  for  measuring 
electrostrictive  coefficients  after  Bohatv 
and  Hauss'iihl , 1977  . 
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Clearly,  If  this  null  condition  can  be  precisely  established,  from  the 
AC  voltage  applied  to  the  standard  piezoid  (d)  and  its  known  dimensions,  the 
amplitude  of  the  strain  at  c can  be  calculated.  Since  this  must  at  null  be 
exactly  the  same  as  the  160  Hz  strain  at  (b)  and  the  voltage  applied  to  (a) 
is  known,  the  appropriate  electrostriction  constant  for  (a)  can  now  be  calcu- 
lated. To  sense  the  null  condition  in  the  capacitor  be,  the  capacitor  is  used 
as  the  tuning  element  of  a 100  MHz  oscillator,  detecting  the  modulation  of  be 
by  a discriminator  and  phase  locked  detector  locked  to  the  160  Hz  driving  oscil- 
lator . 

In  the  system  which  we  have  designed  and  constructed,  the  mode  of  sensing 
the  capacitance  is  different,  and  is,  we  believe,  a radical  improvement  over 
the  Band H design.  For  a simple  parallel  plate  capacitor  dilatometer  operating 
in  dry  nitrogen  the  sensing  capacitance  C is  given  by 


where  £q  is  the  permittivity  of  free  space,  £ the  permittivity  of  dry  nitrogen 
at  atmospheric  pressure,  A the  plate  area,  and  Z the  separation.  The  capacitance 
change  AC  associated  with  a small  change  in  Z (AZ)  will  be  given  by 

£ £A 

AC  = (-  ~-)AZ  = KAZ 

J r 

Obviously  the  sensitivity  constant  K can  be  made  very  large  if  Z can  be 
made  very  small.  In  a practical  system,  however,  a limiting  value  of  Z is  set 
by  slow  dimension  changes  due  to  thermal  drift,  and  capacitance  changes  asso- 
ciated with  gas  pressure  changes,  which  can  take  C^  outside  the  linear  range 
of  any  very  high  sensitivity  discriminator.  In  the  AC  dilatometer,  however, 
since  only  the  component  of  the  plate  separation  AH(2w)  at  the  piezoelectric 


drive  frequency  is  of  interest,  a DC  servo,  actuated  by  an  error  signal  from 
the  static  capacitance  value  can  be  used  to  stabilize  C . 
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I 

A simple  block  diagram  ot  the  full  measuring  system  is  given  in  Fig.  2.2. 

The  components  of  the  mechanical  circuit  inside  the  dotted  box  (a,b,c,d)  are 
identical  in  principle  to  those  used  by  B and  H.  However,  to  sense  and  control 
the  capacitor  be  a General  Radio  1620  capacitance  measuring  assembly  is  used 
(j,k,l). 

To  obtain  maximum  sensitivity,  the  bridge  is  operated  at  a frequency  of 
1 KHz, the  standard  piezoelectric  crystal  is  driven  at  7.5  Hz  and  the  electro- 
strictor  under  study  at  3.75  Hz.  In  the  1620  assembly  the  type  1238  detector 
(2.)  provides  synchronous  detection  of  both  the  in-phase  and  the  quadrature  com- 
ponents of  the  bridge  unbalance  signal,  phase  locked  to  the  1316  bridge  driving 
oscillator  (K) . The  in-phase  DC  signal  output  is  used  to  drive  a Burleigh 
RC42  Servo  amplifier  which  actuates  three  parallel  piezoceramic  mechanical 
actuators  controlling  the  static  plate  separation  (be) . 

A second  phase  locked  detector  (n)  (Keithley  Model  822  with  821  phase 
shifter)  also  looks  at  the  in-phase  unbalance  signal  but  now  with  a time  cons- 
tant X <<  0.1  seconds.  The  output  of  this  detector  is  passed  to  a Princeton 
HR8  lock-in  amplifier  which  is  phase  locked  to  the  7.5  Hz  slave  oscillator 
driving  the  standard  piezoelectric  crystal. 

The  DC  output  of  this  lock-in  provides  the  null  detector  which  senses 
the  condition  under  which  the  plates  b and  c are  being  driven  in  phase  at  the 
same  amplitude,  i.e.  that  the  separation  be  contains  no  component  in  phase 
with  the  driving  oscillators  (a)  and  (d)  . 

2.2.2  Construction  of  the  Dilatometer 

The  mechanical  system  is  built  around  a Burleigh  Model  RC  110  1R  Fabry- 


Perot  Interferometer.  The  IR  interferometer  already  incorporates  PZT  piezo- 
electric drivers  which  actuate  the  movable  interferometer  plate  mount  to  pro- 
vide 8 pmeters  of  linear  movement  and  3.5  pm  alignment  tilt.  This  mechanical 


Block  Diagram 
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frame  (Fig. 2.3)  is  constructed  from  super  invar  (a  0.36  x 10  h/°C)  and 
provides  a 0 - 150  mm  cavity  into  which  to  mount  the  two  driver  crystals  and 
the  sensing  capacitor  with  all  tilt  and  alignment  adjustments  already  built 
in . 

The  elaton  frame  fits  within  a model  RC  75  temperature  control  chamber 
which  stabilizes  the  structure  to  < 0.05°C. 

The  arrangement  of  the  driver  generators  and  sensing  capacitor  is  shown 
schematically  in  Fig.  2.4.  The  standard  piezoelectric  crystal  is  a large  X 
cut  quartz  plate  mounted  between  lucite  supports  (Section  a). 

The  sensing  capacitor  plates  are  of  fused  silica,  with  sputtered  gold 
surfaces  (b)  and  the  electrostrict ion  section  c which  contains  the  electro- 
strictive  plate  under  study  is  a similar  sandwich  structure.  In  this  section 
the  perspex  insulator  on  the  live  side  of  the  electrostr ic tor  crystal  is  made 
very  thick  so  that  the  field  across  this  section  is  low  and  the  electrostric- 
tion  of  the  perspex  does  not  contribute  to  the  total  linear  deformation.  Each 
section  of  the  assembly  is  separately  electrostatically  screened,  and  all  con- 
nections are  through  coaxial  cables. 

Both  the  standard  quartz  driver  and  the  electrostriction  crystal  are  held 
in  place  by  a very  thin  film  of  Eastman  910  cyanoacrylate  ester  cement. 

The  Burleigh  RC  75  thermal  enclosure  which  houses  all  mechanical  components 
is  mounted  on  a EP0I  vibration  isolation  platform,  which  is  itself  supported 
on  an  air  mounted  granite  slab  (Fig.  2.5). 

2.2.3  Dilatometer  Performance 

One  advantageous  feature  of  the  present  design  of  AC  dilatometer  is  that 
the  performance  can  be  verified  stage  by  stage.  First  the  drift  correction  of 


ij 


the  DC  Servo  can  be  checked . 


Fig.  2.3  Burleigh  Model  RC  110  I.R.  Fabry  Perot  cavity 
with  AC  dilatometer  sections  in  place. 

(a)  Standard  piezoelectric  driver. 

(b)  Shielded  dilatometer  capacitor. 

(c)  Electrostriction  driver  section. 


L 


iple  Holder 


Fig.  2.4  Schematic  of  the  dilatometer  and  support  structure 


Fig.  2.5  Mounting  arrangement  for  the  AC  dilatometer 
housed  in  the  Burleigh  RC  75  temperature 
controlled  enclosure. 


After  the  system  has  stabilized  for  some  4 hours,  the  DC  Servo  system  is 
disabled  and  the  intrinsic  drift  rate  of  the  capacitance  determined  by  measuring 
the  bridge  unbalance  signal  as  a function  of  Lime  (Fig.  2.6).  By  rebalancing 
the  bridge  periodically  the  capacitance  increment  for  a given  unbalance  voltage 
can  be  determined,  and  from  the  AC  a corresponding  Ail  can  be  calculated. 

It  is  evident  from  Fig.  2.6  that  without  compensation  the  drift  rate  corres- 
ponds to  a plate  separation  change  of  lX/minute.  Maintaining  the  same  conditions, 
the  DC  Servo  is  now  turned  on,  and  at  the  same  gain  setting  no  drift  is  discern- 
ible at  all.  Turning  the  gain  up  by  a factor  of  100,  however  (Fig.  2.6  right- 

hand  scale),  it  is  evident  that  the  drift  rate  has  been  decreased  by  a factor 

-3  ° 

^2  x 10  and  now  corresponds  to  a linear  dimension  change  '±'.002A/minute . 

Thus  the  detector  amplifier  can  be  maintained  at  110  dB  over  long  periods 
of  time  without  drifting  out  of  the  linear  range  of  the  amplifier. 

To  check  the  AC  sensitivity  the  standard  quartz  crystal  was  driven 
by  a Wavetek  Model  142  oscillator  at  7.45  Hz.  The  conditions  recorded  in 
Table  2.1  were  set  up. 

With  20  volts  RMS  applied  to  the  bridge,  the  lock-in  output  as  a function 

of  voltage  applied  to  the  quartz  crystal  is  given  in  Fig.  2.7.  The  linear 

relation  expected  for  a piezoelectric  is  clearly  evident.  Taking  a value  of 

d = +2.3.10-10  and  a resolution  limit  of  ±0.2  mV  gives  for  the  smallest  resolv- 

-3° 

able  increment  AJ,  = ±5.10  A. 

As  a further  check  of  the  system,  the  lock-in  output  for  a fixed  value  of 
drive  voltage  to  the  quartz  (10.8  V RMS)  was  recorded  as  a function  of  voltage 
applied  to  the  bridge.  If  the  measured  signal  comes  only  from  the  bridge  un- 
balance, this  should  scale  linearly  with  the  applied  voltage  as  is  evident  in 
Fig.  2.8. 

Experiments  are  now  in  progress  to  test  the  electrostrictive  cell  which 


is  the  final  component  in  the  system. 


(GAP) 
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Fig.  2.6  DC  plate  separation  change  AZ  as  a function  of  time  with  and  without 
the  drift  compensation  Servo  amplifier. 
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Table  2.1  Conditions  used  in  the  determination  of  AC  sensitivity 


1.  Temperature 

2.  Capacitance  Bridge: 

Gain 

Capacitance 

Conductance 

3.  Quartz  Driving: 

4.  Band  Pass  Filter: 

5.  Lock-in  Amplifier: 

Signal  Channel 
Mode 

Q 

Time  Const. 

Reference  Channel 
PSD 


29.4  ± 0.04°C 

1 kHz,  20-30  Vrms 
80-100  dB 

250  pF  (plate  gap  70  pm) 

1.4  x 10~4  p V 
7.45  Hz,  0-10  Vrms 

2.5  - 13  Hz 

7.45  Hz 

Band  Pass 
100 

10  sec 

Normal 


k 


\ 
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Deformation  Arrmt  (A) 

0 0.05  0.1  0.15  0.2  0.25 


Fig.  2.7  Dilatometer  output  as  a function  of  driving  voltage  applied  to 
the  quartz  standard  crystal  and  of  equivalent  AC  displacement 
amp  1 i Lude . 


Lock-in  Amplifier  Output  (mV) 


25 


20 


15 


Bridge  Oscillator  Voltage  (Vrms) 


Fig.  2.8  Dilatometer  output  as  a function  of  bridge  oscillator  voltage 
for  fixed  10.8  V rms  applied  to  the  quartz  crystal. 
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2 . 3 Lattice  Theory  of  Elec t rust r i c t i o n 

2.3.1  Technical  Problem 

As  shown  phenomenologically  the  piezoelectric  constants  of  ferroelectr ics 
derived  from  centric  prototypes  (the  case  for  all  present  piezoceramics)  are 
proportional  to  the  electrostriction  constants  and  the  spontaneous  polarization. 
Thus  for  systematic  attempts  to  find  new  materials  with  larger  electromechanical 
coupling  factors  the  microscopic  parameters  which  determine  these  properties 
must  be  well  understood.  A critical  analysis  of  previous  theoretical  treatments 
of  electrostriction  shows  that  all  are  inadequate  even  for  such  simple  solids 
as  alkali  halides.  The  purpose  of  the  present  work  is  to  develop  a more  adequate 
theory  of  electrostriction,  to  apply  it  to  alkali  halides  and,  if  successful, 
to  selected  piezoelectric  and  ferroelectric  materials  occurring  in  simple  struc- 
tures. It  is  hoped  that  the  results  will  provide  heuristic  guidelines  for  maxim- 
izing the  electrostriction  constants  for  a given  crystal  structure. 

2.3.2  Methodological  Approach 

All  previous  theories  of  electrostriction  are  based  either  on  the  rigid 
ion  model,  or  on  various  versions  of  the  shell  model  of  lattice  dynamics.  Obvi- 
ously, the  rigid  ion  model  is  inadequate  because  electronic  polarizability  is 
neglected  altogether.  In  the  available  shell  model  treatments  electronic  polari- 
zability is  only  included  in  the  harmonic  approximation,  in  which  the  strain 
dependence  of  the  electronic  polarizability  is  neglected.  In  addition,  other 
simplifying  assumptions  of  dubious  validity  are  made. 

The  present  work  consists  of  two  phases.  In  the  first  phase  several 
anharmonic  shell  models  are  being  investigated  in  a rigorous  and  consistent 
manner.  Promising  models  which  can  satisfactorily  account  for  electrostriction 
in  alkali  halides  will  be  subsequently  applied  to  piezoelectric  and  ferro- 
electric materials  occurring  in  simple  crystal  structures.  For  the  second 
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phase  empirical  pseudopotential  calculations  in  conjunction  with  the 
dielectric  screening  approach  are  envisaged  which  would  relate  the  empirical 
parameters  in  the  anharmonic  shell  models  to  more  basic  microscopic  physical 
quantities . 

2.3,3  Technical  Results 

2. 3. 3.1  Definition  of  Electrostriction  Constants 

Of  the  various  free  energy  functions  pertaining  to  the  elastic 


dielectric  we  choose  the  Helmholtz  free  energy  F(T,n..,E  ) as  a function  of 

tj  k 

temperature  T,  strain  n . . , and  electric  field  E,  (i, j ,k=l ,2 , 3) . The  electro- 

ij  k 

strictive  constants  are  defined  as  the  third  derivative  of  the  free  energy 
according  to 


2f  = ( XX \ 

ijkl  v3n  SE^E^T 


with  the  definitions  for  the  stress  tensor 


, 3F  . 

ij  - (3nlj)T,Ek 


and  of  the  dielectric  constant  tensor 


(2.1) 


(2.2a) 


( ii- . > 

V 3E.  3E/T,n.  . 

k l ij 


(2.2b) 


one  obtains  from  equ  (2.1)  after  reversing  the  sequence  of  the  derivatives  the 
alternative  expression 


3e 


3 -‘-M., . - - <--■) 
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3n-  . T,E 

lj  m 


(2.3) 


The  electrostrictive  constants  have  the  symmetry  properties 
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According  to  eau  (2.3)  the  electrostrictive  constants  i ...  , are  the  first  order 

ljkl 

strain  derivatives  of  the  dielectric  constant  tensor.  From  equ  (2.1)  it  is 
apparent  that  the  electrostrictive  constants  depend  on  the  third  order  coupling 
parameters  among  the  atoms  and  the  electronic  degrees  of  freedom  and  therefore 
represent  first  order  anharmonic  effects.  Therefore  lattice  theoretical  models 
of  electrostrict ion  must  include  anharmonic  effects  to  first  order  in  a consist- 
ent manner. 

£.3.3.2  Shell  Models  in  Lattice  Dynamics 

The  shell  model  was  introduced  almost  two  decades  ago  (Cochran,  1959; 
Woods  et  al,  1960;  Cochran,  1961;  Dick  and  Overhauser,  1958;  Hardy,  1962)  and 
is  now  the  most  widely  used  theoretical  model  in  the  lattice  theory  of  ionic 
and  covalent  crystals  (See,  for  example,  the  review  articles  by  Cochran,  1971, 
Bilz  et  al.,  1974,  and  Hardy,  1974).  The  reasons  for  its  success  are  that  it 
provides  a simple,  but  sufficiently  adequate  method  for  treating  effects  asso- 
ciated with  the  electronic  polarizability  of  the  ions,  and  that  it  explicitly 
takes  into  account  the  difference  in  the  location  of  the  electron  shells, 
which  determine  the  interatomic  short-range  exchange  and  overlap  forces,  and 
the  location  of  the  nuclei,  which  are  associated  with  the  inertial  forces. 

This  is  achieved  by  approximating  the  electronic  charge  distribution  by  (rigid) 
spherical  shells  which  are  coupled  by  means  of  springs  to  the  ion  cores  composed 
of  the  nuclei  and  the  inner  electron  shells. 

In  the  rigid  ion  model  (Kellermann,  1940)  the  interatomic  forces  are  com- 
pletely specified  by  the  ionic  charges  Z e (e  = electronic  charge;  < labels 
the  different  types  of  ions),  which  give  rise  to  the  long-range  Coulomb  forces, 
and  the  derivatives  of  the  short-range  interaction  potential  i(R).  The  latter 
are  conveniently  described  by  the  dimensionless  parameters: 
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(2.5a) 


(2.5b) 
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(2.5c) 


R^  6.=  1,2,...)  denotes  the  interatomic  distances  between  first,  second,  etc. 
neighbors,  and  V the  volume  per  unit  cell.  Alternatively,  for  the  Born-Mayer 
form  of  the  interatomic  potential, 

4>  (R)  = be-R/p  (2.6) 

the  two  parameters  b and  p are  required.  One  of  the  parameters  in  (2.6),  or 
one  of  the  parameters  B^  according  to  equ  (2.5a)  may  be  eliminated  by  means  of 
the  equilibrium  condition.  Noncentral  forces  may  be  described  if  the  two 
parameters  A^  and  B^  are  taken  as  independent,  so  that  they  describe  the  second 
derivatives  of  the  interatomic  potential  for  displacements  parallel  and  per- 
pendicular to  the  direction  connecting  two  ions  according  to 

A=  (^)(4| 
e dR  1 

B = (2V)  (^|)  , 
e 9rZ  -L 


(2.7a) 

(2.7b) 


In  the  original  and  simplest  version  of  the  shell  model,  the  rigid  shell 
model,  two  additional  parameters  are  required  for  each  polarizable  ion.  They 
are  the  shell  charges  and  the  core-shell  force  constants  which  enter 
the  expressions  for  the  electronic  polarizability.  In  later  versions  of  the 
shell  model  additional  electronic  degrees  of  freedom  are  introduced,  or  non- 


linear  effects  are  included.  For  example,  in  the  breathing  shell  model 
(Schroder,  1966;  Niisslein  and  Schroder,  1967)  a radial  deformation  of  the  electron 
shell  is  allowed.  In  the  charge  transfer  model  (Verma  and  Singh,  1969;  Verma 
and  Agarwal,  1973;  Feldkamp,  1972)  charge  transfer  between  ions  is  introduced, 
which  gives  rise  to  three-body  forces.  In  the  deformable  ion  model  (Jaswal 
and  Dilly,  1977;  Jaswal,  1978)  the  deformation  polarizability  and  the  electrical 
polarizability  are  taken  as  different.  Among  the  major  achievements  of  these 
more  elaborate  shell  models  are  that  they  provide  better  agreement  between 
calculated  and  measured  phonon  dispersion  curves  along  different  crystallographic 
symmetry  directions  with  fewer  empirical  parameters  than  the  rigid  shell  model  , 
because  the  short-range  interactions  need  to  be  considered  only  up  to  fewer 
neighbors.  For  example,  for  MgO  the  charge  transfer  model  gives  excellent 
agreement  with  six  parameters  (Verma  and  Agarwal,  1973),  and  for  ZnS  (zinc- 
blende)  only  five  parameters  are  required  (Jaswal,  1978). 

In  the  charge  transfer  model  (Verma  and  Singh,  1969;  Verma  and  Agarwal, 

1973;  Feldkamp,  1972)  charge  transfer  between  nearest  neighbor  ions  is  intro- 
duced as  a function  of  the  interionic  distance  and  is  described  by  a function 
f(R).  The  additional  parameters  reauired  in  the  theory  are  f = f(R^); 

f'  = (df/dR)„  and  f"  = (d2f/dR2)„  . 

K K 

O O 

All  the  parameters  of  the  shell  model  discussed  above  are  required  for 
the  description  of  properties  in  the  harmonic  approximation,  with  the  exception 
of  CK  as  defined  in  equ  (2.5c)  and  the  charge  transfer  parameter  f"  mentioned 
in  the  previous  paragraph,  which  are  required  for  the  description  of  first 
order  anharmonic  effects,  such  as  the  electrostrictive  constants.  The  constants 
CK  may  be  either  treated  as  independent  empirical  parameters,  or  calculated 
from  the  parameters  A^,  IK  by  assuming  the  Born-Mayer  form  (2.6)  for  the 


interatomic  potential. 
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2.3.3. 3 Previous  Theories  of  Electrostr iction 

In  Table  2.2  the  theoretical  models  are  specified  which  have  previously 

been  used  to  account  for  the  electrostrictive  constants  of  alkali  halides. 

In  the  rigid  ion  model  of  Grindlay  and  Wong  (1969),  denoted  by  RI-1,  short- 

range  interaction  is  considered  between  nearest  neighbors  only,  and  the  Born- 

Mayer  potential  (2.6)  is  assumed.  The  only  free  parameter  of  this  model  is 

the  hardness  parameter  p which  is  fitted  to  the  static  dielectric  constant  t . 

o 

For  comparison  a second  rigid  ion  model,  denoted  by  RI-2,  is  included  in  Table  2.1 
in  which  the  first  neighbor  interaction  is  to  be  of  noncentral  type  up  to  second 
order  in  the  ionic  displacements,  in  which  second  neighbor  interactions  of 
central  force  type  are  included,  and  in  which  the  form  of  the  interatomic  poten- 
tial is  not  specified.  The  three  harmonic  parameters  A^,  B^,  k^  are  fitted  to 

the  three  elastic  constants,  and  the  two  anharmonic  parameters  are  fitted  to 

I 

the  first  pressure  derivatives  of  two  of  the  elastic  constants. 

The  shell  model  (SM)  included  in  Table  2.2  is  a generalized  version  of 
the  model  considered  by  Srinivasan  and  Srinivasan  (1972).  It  differs  from  RI-2 
in  that  the  polarizability  of  the  anions,  described  in  terms  of  the  two  additional 
parameters  Y and  k,  is  included,  and  in  that  the  five  harmonic  parameters  are 
determined  from  the  three  elastic  constants,  from  the  static  and  optical  diel- 
ectric constants  and  respectively,  and  from  the  transverse  optical  frequency 

u)T0' 

In  the  model  of  Goyal  et  al.  (1978),  denoted  by  SM-MB , three-body  forces 
arising  from  the  charge  transfer  mechanism  are  included,  but  short-range  forces 
are  considered  only  between  first  nearest  neighbors,  are  assumed  to  be  of 
central  force  type,  and  are  described  by  a Born-Mayer  potential.  Therefore, 
the  only  anharmonic  parameter  is  the  second  derivative  f",  which  is  fitted  to 
the  pressure  coefficient  of  the  bulk  modulus. 


Table  2.2  Previous  theoretical  models  used  for  calculation  of  electrostrive  constants 
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Rigid  ion  model;  SM  = shell  model;  MB  = many  body  force  model 
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In  Table  2.3  the  experimental  el ectrostr ietlve  constants  in  Voigt 
notation*  and  the  pressure  coefficient  of  the  static  dielectric  constant** 
are  compared  with  the  results  calculated  from  the  four  models  for  two  selected 
alkali  halides.  In  order  to  eliminate  variations  arising  from  different  input 
parameters  the  values  recalculated  by  us  from  identical  input  parameters  are 
given  instead  of  the  values  calculated  by  the  original  authors.  The  experi- 
mental values  are  from  Bohaty  and  Haussuhl  (1978)  . 

It  is  apparent  that  for  NaCl  model  SM-MB  gives  the  best  overall  agreement 
and  the  correct  sign  for  f . However,  for  f^  model  RI-2  gives  the  best 
agreement,  and  the  value  of  f^  for  model  SM-MB  is  far  too  small.  For  KBr  model 
RI-1  gives  the  best  overall  agreement.  However,  only  models  SM  and  SM-MB  give 
the  correct  sign  of  f,^,  although  again  the  calculated  coefficients  are  far 
too  small.  For  NaCl  model  RI-2  gives  the  best  agreement  for  f^*  and  for  KBr 
both  models  RI-1  and  RI-2  give  the  best,  but  still  very  poor  agreement  for  f ^ • 

The  relatively  good  agreement  (which  is,  though,  still  far  from  satisfactory) 
obtained  for  KBr  with  the  rigid  ion  model  RI-1  is  surprising  and  must  be  attri- 
buted to  a cancellation  of  errors,  because  the  model  is  the  simplest  of  all 
models  considered  and  does  not  include  electronic  polarizability.  Therefore 
the  optical  dielectric  constant  would  be  equal  to  one,  and  all  piezooptic 
constants  would  be  zero. 


*Index  pairs  (ij)  occurring  in  tensor  notation  are  replaced  by  single  indices 
p according  to  the  scheme  (ij)  = 11,22,33,23  and  32,13  and  31,  12  and  21 


M = 1»  2 , 3,4, 5, 6 . 

**This  quantity  is  related  to  the  electrostrictive  constants  defined  by  (2.1)  and 
(2.3)  according  to 

<-l|)  ‘ <3f)(fn  + 2fi2> 


where  B denotes  the  bulk  modulus. 


Model 

fll 

f 12 

f 44 

Experiment 

-2.67 

0.62 

-0.97 

RI-1 

-3.03 

0.36 

0.17 

RI-2 

-5.03 

0.55 

0.35 

SM 

-2.23 

0.21 

0.009 

SM-MB 

-2.56 

0.19 

-0.03 

Experiment 

-2.56 

0.82 

-0.66 

RI-1 

-2.57 

0.31 

0.15 

RI-2 

-5.32 

0.31 

0.15 

SM 



-1.98 

0.08 

-0.09 

SM-MB 


-2. 47 


0.11 


-0.08 
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Model  SM-MB  is  the  most  elaborate  and  realistic  of  the  models  considered 
in  that  electronic  polarizability  effects  are  included  and  even  some  form  of 
anharmonicity  of  the  electron  shells  is  taken  into  account  through  the  charge 
transfer  parameter  f".  However,  as  will  be  discussed  in  Section  2. 3.3.4,  the 
anharmonic  shell  effects  are  treated  in  an  inconsistent  manner  in  SM-MB.  In 
addition,  other  potentially  important  effects  are  omitted  in  SM-MB,  such  as 
nonlinearity  of  the  core-shell  coupling,  and  second  neighbor  short-range  repulsive 
forces.  Perhaps  the  relatively  poor  agreement  of  this  model  for  KBr  could 
also  be  attributed  in  part  to  the  neglect  of  cation  polarizability. 

Clearly,  theoretical  work  on  electrostr iction  in  piezoelectric  and  ferro- 
electric crystals  can  only  be  started  after  an  adequate  model  and  an  adequate 
theory  for  electrostriction  in  alkali  halides  are  available.  The  objective 
of  the  work  pursued  under  the  present  contract  is  to  develop  such  a theory  and 
to  provide  this  theoretical  framework  for  subsequent  application  to  piezo- 
electric and  ferroelectric  materials. 

2. 3. 3. 4 Lattice  Dynamics  of  Anharmonic  Shell  Models 

The  lattice  dynamical  framework  for  the  calculation  of  numerous  anhar- 
monic properties  has  been  derived  and/or  reviewed  by  Cowley  (1964,1965). 

However,  as  first  pointed  out  by  Oitmaa  (1967)  the  anharmonicity  included  by 
Cowley  and  most  of  the  subsequent  applications  refers  only  to  the  nonlinearity 
of  the  interionic  forces  with  respect  to  ionic  displacements,  but  the  anhar- 
monicity associated  with  the  intraionic  forces  is  neglected,  that  is  the  core- 
shell coupling  is  harmonic.  Notable  exceptions  are  calculations  of  thermal 
expansion  based  on  the  quasiharmonic  approximation  (See  for  example  the  review 
article  by  Clyde  and  Klein,  1971),  of  the  second-order  Raman  spectra  of  alkali- 
earth  oxides  (Rieder  et  al.,  1973;  Migoni  et  al . , 1976),  and  of  the  mode; 
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softening  in  KTaO.^  (Migoni  et  al.,  1978).  The  thermal  expansion  calculations 
are  based  on  a volume  dependent  core-shell  force  constant,  and  in  the  theory 
of  the  second  order  Raman  spectra  and  of  the  mode  softening  a quartic  term  in 
the  displacements  is  introduced  in  the  energy  for  the  core-shell  displacements 
and  is  found  to  be  of  crucial  importance. 

Since  the  electrostrictive  coefficients  describe  the  strain  dependence  of 
the  static  dielectric  constant  tensor,  and  since  the  dielectric  constant 
includes  electronic  polarizability  effects  it  would  appear  mandatory  to  include 
intraionic  anharmonicity  in  any  adequate  theory  of  electrostriction . As  men- 
tioned in  Section  2.3.3. 2 the  only  model  in  which  this  is  attempted,  SM-MB,it  is 
done,  however,  in  a manner  which  is  inadequate  for  the  following  three  reaons:  j 

First,  the  primary  feature  of  intraionic  anharmonicity,  the  nonlinear  core-shell  j 

coupling  is  not  taken  into  account.  Second,  the  anharmonicity  resulting  from 
the  charge  transfer  mechanism  primarily  affects  the  interionic  forces,  but  not 
the  intraionic  forces.  Third,  the  anharmonic  effects  are  introduced  in  the 
equations  derived  by  Srinivasan  (1968) , which  refer  to  a shell  model  with  purely 
harmonic  intraionic  forces.  As  pointed  out  by  Oitmaa  (1967),  the  effective 
anharmonic  coupling  parameters  for  a fully  anharmonic  shell  model  depend  in  an 
essential  manner  on  the  intraionic  anharmonicity. 

In  the  present  work  the  effective  coupling  parameters  have  been  derived 
consistently  for  a fully  anharmonic  shell  model,  including  both  inter  ionic  and 
intraionic  anharmonicity.  With  these  correct  expressions  exact  and  general 
expressions  for  the  pertinent  first  order  anharmonic  coefficients  are  being 
derived  at  present.  This  includes  the  third  order  elastic  constants,  the  photo- 
elastic coefficients,  and  the  electrostrictive  coefficients.  Application  to 
simple  crystal  structures,  especially  the  rocksalt,  zincblende  and  perovskite 
structures,  is  planned  for  the  future. 
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2.3.4  Summary  and  Conclusions 

All  previous  theoretical  treatments  of  the  electrostrictive  constants 
are  inadequate  from  a theoretical  point  of  view  and  because  of  poor  agreement 
with  the  experimental  data  for  alkali  halides.  We  are  in  the  process  of 
deriving  consistent  and  general  equations  for  the  first  order  anharmonic  mate 
rial  constants,  including  the  electrostrictive  constants,  in  which  anharmonic 
effects  associated  both  with  interionic  and  intraionic  forces  are  included. 
Considerable  improvement  of  the  agreement  between  theory  and  experiment  is 
expected  for  alkali  halides,  and  subsequent  application  to  piezoelectric  and 
ferroelectric  materials  is  planned. 
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3.  PIEZOELECTRIC  COMPOSITES 
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3 . 1 Introduction 

Earlier  studies  on  ONR  contract  N000.1 4-76-C-0515  have  demonstrated  quite 
clearly  that  for  voltage  generating  devices  such  as  hydrophones,  where  the 
figure  of  merit  for  the  piezoelectric  material  involves  the  hydrostatic  voltage 
coefficient  g^,  or  a produce  g^d^,  composite  piezoelectrics  can  be  fabricated 
which  have  properties  which  are  superior  to  those  of  either  single  component. 

In  the  earlier  study,  a ceramic : plast ic  composite  of  controlled  phase 
interconnection^'*^  was  fabricated  by  a ceramic  replication  process,  using 

(3.2) 

the  three  dimensionally  porous  coral  porites: porites  as  a template  structure 
The  reasons  for  seeking  to  control  the  phase  interconnection  in  diphasic  systems 
which  try  to  exploit  coupled  properties,  as  in  the  electromechanical  transducer, 
have  been  given  in  detail  in  3.1.  The  replication  process  by  which  coral  was 
exploited  as  a template  microstructure  to  produce  the  required  form  of  phase 
interconnection  for  a transducer  with  high  hydrostatic  sensitivity  has  been 
given  in  3.2. 

A significant  disadvantage  of  the  replamine  process  described  in  3.2  as 

a mode  for  reproducing  a 3:1  connected  structure,  is  that  it  is  topologically 

impossible  to  use  the  template  more  than  once,  i.e.  a new  block  of  the  template 

material  of  suitable  size  and  shape  is  required  in  the  fabrication  of  each 

(3.3) 

device.  Both  theoretical  studies  in  this  Laboratory  and  recent  macro- 

(3.4) 

modeling  studies  at  NRL  ' have  confirmed,  however,  that  the  ideal  phase 
connection  for  the  type  of  transduction  required  is  really  3:1  and  not  full  3:3. 

In  the  replamine  studies  this  was  also  apparent,  as  the  coral  porites: 
goniopera  which  has  a preferred  anisotropy  in  its  structure  with  strong  connec- 
tions in  1 dimension,  and  weaker  links  in  the  two  orthagonal  directions  gave 
better  results  than  the  more  isotropic  porites: porites . 


Fortunately,  3:1  connection  is  in  many  ways  more  simple  to  architect  in 
a composite  than  is  the  3:3  structure  of  the  coral.  In  these  new  studies, 
two  methods  are  being  explored  for  the  fabrication  of  ceramic : plastic  composites 
with  dominant  3:1  connectivity  and  a range  of  scales  in  the  microstructure. 

The  first  is  a modification  of  the  original  replamineform  process,  but 
using  now  a manufactured  prototype  structure  of  fine  steel  rods  on  regularly 
spaced  centers.  Again  a wax  negative  is  made  from  the  master  template,  but 
now  after  the  wax  has  set  the  template  can  be  extracted  and  reused  many  times. 

In  the  second  effort,  which  started  over  the  last  three  months,  very  fine 
PZT  rod  is  being  extruded.  These  PZT  fibers  appear  to  offer  intriguing  possi- 
bilities for  manipulation  in  the  flexible  green  state  to  form  a wide  range  of 
composites  with  3:1  and  with  other  more  complex  modes  of  interconnection. 

3 . 2 Casting  from  Metal  Template  Structures 

The  basic  idea  was  to  generate  a structure  of  regularly  spaced  columns  of 
PZT  supported  on  a PZT  base  plate,  on  a scale  where  the  columns  could  be  down 
to  100  pmeters  in  diameter  and  arranged  in  a regular  two-dimensional  pattern. 

The  most  obvious  template  for  this  type  of  structure  is  a wire  brush,  with 
bristles  of  the  required  diameter  arranged  on  regularly  spaced  centers.  After 
the  wires  have  been  coated  with  a suitable  parting  agent,  wax  can  be  cast  around 
the  wires.  When  the  wax  has  set,  the  wires  can  be  withdrawn,  leaving  a wax 
template  with  regularly  spaced  columnar  holes.  The  holes  can  now  be  filled 
with  a PZT  slip  by  vacuum  impregnation,  the  slip  dried,  the  wax  melted  out  and 
the  ceramic  fired,  just  as  in  the  replamine  process. 

The  major  advantage  would  be  that  the  template  is  not  destroyed,  but  is 
removed  intact,  so  that  an  expensive  carefully  prepared  template  could  be 
reused  many  times  in  the  process.  Since  the  scale  of  wire  used,  and  the  spacing. 
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are  open  to  choice,  intentionally  periodic  structures  could  be  fabricated 
with  the  scale  and  the  periodicity  controllable  in  the  plane  of  the  template. 

The  flow  sheet  for  the  processing,  which  is  closely  similar  to  that  of 
the  replamine  processing  is  given  in  Fig.  3.1. 

Initial  experiments  were  carried  out  using  a low-density  structure  con- 
sisting of  8 steel  rods  300  pmeters  in  diameter  1 cm  long,  spaced  around  the 
perimeter  of  a 5 mm  diameter  circle.  The  casting  wax  used  was  a Kerr  inlay 
wax,  though  subsequent  experiments  have  used  a number  of  different  natural  and 
synthetic  waxes.  Investment  of  the  wax  megatives  was  with  a PZT  501  slip, 
and  the  integrity  of  the  invested  PZT  columns  was  checked  by  x-ray  shadow 
photography. 

A major  problem  with  this  method  occurs  at  the  stage  of  wax  removal.  The 
meltdown  operation  is  exceedingly  critical,  the  melting  wax  breaking  up  the 
fragile  PZT  columns.  Once  the  wax  is  removed  there  appeared  to  be  no  serious 
problem  in  firing  the  PZT  columns  without  serious  warpage. 

Experiments  are  now  proceeding  along  the  following  lines: 

1.  Improving  the  drying  of  the  PZT  slip  by  using  vacuum  drying  and 
microwave  heating. 

2.  Strengthening  the  PZT  columns  by  decreasing  the  height  and  thus 
increasing  the  diameter-to-height  ratio. 

3.  Experimenting  with  "abrupt"  melting  synthetic  waxes  which  may  be 
melted  by  radiant  heat  from  the  surface,  i.e.  melting  back  layer  by  layer. 

4.  Testing  campher  based  waxes  which  sublime  directly  without  melting. 

5.  Using  a thin  layer  of  slip  on  both  surfaces  of  the  wax  replica 
(using  a much  higher  density  of  rods)  so  that  the  rods  are  supported  at  both 
ends  as  the  wax  melts  out. 
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CAST  ; Fire  at  1280°C 

DENSE  FIBER 

CERAMIC  j- 

STRUCTURE 

FIBERS 

Fig.  3.1  Flow  chart  for  lost  wax  casting  process  of  replication 
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Our  experience  to  date  suggests  that  the  template  structures  can  be 
generated  and  controlled  quite  easily,  and  offer  a very  wide  range  of  possi- 
bility both  in  scale  factor  and  in  make-up  in  the  plane.  If  the  step  of  wax 
removal  can  be  accomplished  satisfactorily,  the  method  should  make  an  inter- 
esting and  more  controllable  alternative  to  the  present  replamine  processing. 

3.3  Extrusion  of  PZT  Fibers 

An  attractive  alternative  to  the  processing  discussed  above  would  clearly 
become  possible  if  the  PZT  itself  could  be  fabricated  as  a fine  fiber.  Pre- 
liminary studies  with  a ratiier  simple  extrusion  press  suggest  that  this  may 
in  fact  be  possible  for  an  extensive  range  of  fiber  diameters. 

The  extrusion  apparatus  which  has  been  used  consists  of  a small  Carver 
Model  M hydraulic  laboratory  press,  a one-inch  diameter  ram  and  die  set  with 
a detachable  end  plate.  The  paste  is  extruded  through  a single  hole  in  the 
end  plate,  and  to  avoid  curling  of  the  extruded  fiber  it  was  found  desirable 
to  use  an  extrusion  channel  which  is  as  long  as  possible.  For  the  present 
experiments  holes  of  40,  30,  20  and  13  thousandths  of  an  inch  (mils)  have  been 
used.  The  40  and  30  mil  holes  were  drilled  in  3/4"  plate,  the  20  and  13  mil 
extrusion  holes  in  5/8"  plate. 

The  paste  used  for  these  studies  was  of  a commercial  PZT  501A  formulation 
and  was  made  up  from  90  wt%  PZT,  8 wt%  water,  and  2 wt%  polyvinyl  alcohol  (PVA) . 
The  paste  was  ball  milled  in  a plastic  lined  mill  for  12  hours  then  extruded 
several  times  through  a 40  mil  die  to  ensure  homogeneity. 

After  the  extruded  fiber  has  dried  in  air,  it  is  cut  into  suitable  lengths, 
and  the  binder  burned  out  over  a 1/2  hour  period  at  550°C.  Sintering  was 
carried  out  at  1310°C  in  a closed  crucible  with  control  of  the  PbO  partial 
pressure.  For  the  single  fibers,  densities  greater  than  95%  of  theoretical 


have  been  obtained.  Both  straight  and  curved  fibers  have  been  fired,  and  it 
appears  that  the  general  shape  imparted  in  the  green  state  can  be  preserved 
through  the  firing  process. 


Work  is  now  in  progress  exploring  various  methods  for  processing  and 
assembling  fiber  based  composites. 
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GRAIN  ORIENTED  CERAMICS 


4 . L Int reduction 

Mixed  bismuth  oxide  compounds  with  layer  type  structures  form  an  exten- 
sive family  of  ferroelectr ics  with  large  spontaneous  polarizations  and  high 
Curie  temperatures.  The  prototypic  symmetry  for  the  family  is  4/mmm,  and 
most  members  have  ferroelectric  phases  in  orthorhombic  mm2  with  only  four 
equivalent  planar  orientations  for  P . While  the  Curie  temperatures,  polari- 
zations, and  electrostrict ion  constants  suggest  that  the  family  should  be  of 
major  interest  for  piezoceramic  applications,  the  large  planar  anisotropy 
prevents  effective  poling  of  a randomly  axed  polycrystal. 

Some  few  years  ago  we  proposed  to  investigate,  in  an  earlier  ONR  program, 
the  possibility  of  generating  grain  oriented  ceramics  with  a dominant  planar 
anisotropy,  so  that  the  poling  field  could  be  always  in  a plane  of  easy 
ferroelectric  switching  and  the  ceramic  could  be  poled.  Unfortunately  the 


resources  were  not  adequate  to  undertake  this  topic. 

(4.1  4 2 4.3) 

In  the  meantime,  recent  studies  in  Japan  ’ ’ 


have  produced 


excellent  grain  oriented  Bi^Ti^O^*  and  other  compositions  in  the  layer  struc- 
ture family  using  a uniaxial  hot-pressing  technique  which  permits  limited 
extrusion.  This  clever  process  produces  most  effective  orientation  in  the 
ceramic  grain  structure  for  samples  fabricated  from  powders  which  exhibit 
little  or  no  shape  anisotropy.  It  is,  however,  an  expensive  process  for 
manufacture. 

In  the  present  program  we  are  exploring  an  alternative  approach  which 
may  be  significantly  less  expensive  and  appears  accessible  for  a wide  range 
of  the  layer  structure  compounds. 


4.2  Fused  Salt  Processing  of  Layer  Structure  Oxides 

Our  approach  is  to  use  a synthesis  step  which  appears  to  permit  the 
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fabrication  of  powders  in  the  bismuth  oxide  layer  structure  family  which 
have  a marked  platey  morphology . These  powders  can  then  be  processed  by 
a range  of  conventional  ceramic  forming  techniques,  sedimentation,  tape 
casting,  slip  casting,  etc.  which  are  known  to  preserve  the  built-in  aniso- 
tropy of  the  powder. 

It  is  an  observation  of  almost  all  attempts  to  grow  single  crystals 
of  the  bismuth  titanate  Bi^Ti^O^  that  whatever  the  flux  system,  the  crystals 
grow  with  a mica-like  habit  and  extremely  platey  morphology.  Thus  we  were 
tempted  to  try  rapid  exsolution  from  a molten  salt  bath  as  a possible  tech- 
nique for  generating  microcrystals. 

The  flux  used  was  a 1:1  mixture  of  KCl:NaCl.  For  Bi,Ti_0,.  stoichiometric 

4 3 12 

proportions  of  the  Bi^O^  and  TiC^  were  mixed  with  the  1:1  salt  in  weight 
ratios  1:1  up  to  1:10.  After  milling  the  mixtures  were  fired  in  a covered 
alumina  or  platinum  crucible  to  a temperature  of  1150°C,  then  cooled  at  the 
free  running  furnace  cooling  rate.  After  cooling,  the  salt  was  dissolved 
out  by  repeated  washing  in  water  and  the  powders  were  dried. 

A typical  scanning  electron  micrograph  of  the  powder  produced  in  this 
manner  is  shown  in  Fig.  4.1,  where  the  very  plate-like  character  is  obvious. 
Contrary  to  our  initial  expectation,  the  average  particle  size  in  the  major 
dimension  did  not  depend  markedly  upon  the  cooling  rate  of  the  flux.  Slow 
cooling  appeared  to  give  rather  more  angular  particles  with  better  developed 
(hkO)  faces,  but  neither  the  major  dimensions  of  the  sheets  nor  the  thick- 
ness was  markedly  increased  by  the  slower  cooling  rates. 

A similar  family  of  experiments  have  been  carried  through  for  the  com- 
pound bismuth  tungstate  (B^WO^.)  . In  this  system  a slightly  lower  maximum 
firing  temperature  appeared  to  produce  optimum  results.  A scanning  electron 
micrograph  of  this  powder  is  shown  in  Fig.  4.2. 
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1 


Fig.  4.1  Bi4Ti30^2  powder  produced 

by  molten  salt  process. 
Scanning  electron  micro- 
scope picture  magnification 
indicated . 
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Fig.  4.2  M2WO6  powder  produced 
by  molten  salt  process. 
Electron  micrograph 
magnification  indicated. 
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Experiments  are  now  proceeding  to  tape  cast  grain  oriented  ceramic  sheets 
from  these  precursor  powders.  The  initial  results  with  Bi^Ti.^O^  look  most 
encouraging  with  x-ray  diffraction  showing  a strongly  preferred  c orientation 
normal  to  the  plane  of  the  fired  sheets. 
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5.  FERROELECTRIC  CRYSTAL  GROWTH 

5 . 1 Introduct Ion 

The  rationale  for  establishing  crystal  growing  capability  on  this  pro- 
gram has  been  discussed  earlier.  Since  the  solid  solution  materials  which  are 
of  most  interest  are  all  necessarily  incongruently  melting,  flux  growth  tech- 
niques are  required,  and  the  initial  focus  has  been  in  this  area.  The  Laboratory 
has  had  a long  history  of  successful  high  temperature  flux  growth  with  other 
oxide  systems,  but  recently  this  activity  has  been  reduced,  so  much  of  the 
work  accomplished  over  the  first  months  of  this  contract  has  of  necessity  been 
concerned  with  refurbishing  and  updating  the  furnaces  and  temperature  control 
equipment  for  these  new  studies. 

In  the  following  a brief  description  is  given  of  the  vertical  de-mountable 
type  furnace  which  is  being  used  so  as  to  permit  the  flux  to  be  decanted  before 
solidification,  and  of  the  new  high  precision  |J  data-track  system  which  has 
been  installed  to  permit  the  precise  program  control  of  two  furnaces.  The  use 
of  this  equipment  to  grow  single  crystals  of  lead  magnesium  niobate  for  electro- 
striction  studies  is  described,  and  the  present  ongoing  program  for  growing 
Pb^Mgl^Og : PbTiO^  and  PbZrO^:PbTiO^  crystals  briefly  reviewed. 

5 . 2 Furnace  and  Control  Equipment 

The  demountable  tube  furnace  is  of  rather  simple  construction.  The  mullite 
inner  sheath,  closed  at  one  end,  is  approximately  21”  long  and  5"  in  diameter. 

This  tube  is  supported  on  a firebrick  base  surrounding  a 3-3/4"  diameter  lower 
entry  hole,  and  is  inside  a rectangular  box  formed  by  high  temperature  insula- 
ting bricks  held  in  a Syndanyo  sheet  chamber.  The  top  of  the  chamber  is  com- 
posed of  removable  bricks  drilled  to  accommodate  four  Crystolon  spiral-cut 


silicon  carbide  heating  elements. 


43 


r 


The  crucible  to  be  heated  is  supported  on  a fixed  alumina  column  which 


enters  through  the  lower  entry  hole.  The  whole  furnace  box  is  mounted  on 
vertical  rails  and  is  counterweighted  so  that  it  lifts  off  the  alumina  pedestal 
allowing  easy  access  to  the  growth  crucible. 

The  temperature  controlling  thermocouple  is  attached  to  the  sheath  tube 
at  the  hot  zone,  and  moves  with  the  furnace.  The  temperature  sensing  couple 
is  at  the  top  of  the  support  pedestal  in  contact  with  the  growth  crucible. 

Power  level  required  to  maintain  the  furnace  at  1300°C  is  about  6KVA 
(30  amps  at  200  volts). 

A stop  on  the  demount  rails  positions  the  crucible  in  the  hot  zone  of  the 
furnace,  but  final  positioning  can  be  trimmed  by  a lab  jack  under  the  pedestal 
structure. 

Temperature  control  for  the  furnace  is  provided  by  a Research,  Inc.  640  U 
temperature  control  unit  regulating  the  output  of  a Payne  80  amp  phaser  unit. 
The  control  program  is  established  by  a new  Research,  Inc.  p data-track  digital 
controller  with  Master  and  Blend  station  set.  These  units, together  with  a 
Fluke  digital  thermometer,  voltmeter  and  ammeter  for  the  furnace  current  and 
the  appropriate  switches,  are  mounted  in  a mobile  rack  unit  (Fig.  5.1). 

The  heart  of  the  control  system  is  the  p data-track  unit.  Working  on  a 
12-bit  digital  coding,  the  unit  provides  a control  voltage  resolution  of 
1/4000.  Using  a Master  and  Blend  station  unit,  however,  this  control  signal 
can  be  arranged  to  cover  anything  from  0 - 100%  of  the  full  temperature  range, 
i.e.  shorter  segments  of  the  temperature  range  can  be  covered  with  very  high 
resolution  as  is  often  required  in  crystal  growth. 

In  the  simplest  form,  the  microprocessor  based  controller  which  is  timed 
by  its  own  internal  clock  breaks  the  program  into  19  linear  segments,  duration 
and  rate  of  change  in  each  segment  being  controllable  over  a very  wide  range. 
With  pre-programmed  ROM  the  resolution  can  be  further  improved  to  127  linear 
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program  segments,  each  with  duration  1 to  1,999  seconds,  minutes,  or  hours. 
Transfer  between  program  segments  can  also  be  used  to  activate  relay  closure, 
and  up  to  seven  external  commands  can  be  generated  by  the  system  for  ancillary 
funi  t ions . 

In  conjunction  with  the  demountable  box  kiln,  the  controller  provides  a 
very  precise  but  highly  flexible  system.  With  the  Master  and  Blend  stations 
two  completely  independent  control  program  signals  can  be  generated,  and  we 
anticipate  adding  a second  furnace  to  complete  the  unit. 

5 . 3 Growth  of  Lead  Magnesium  Niobate 

Lead  magnesium  niobate  Pb (Mg^y 3^*2/ 3^3  crystals  were  grown  from  a flux 
of  excess  PbO  with  B2O.J,  using  a composition  ratio  which  had  been  found  to 
work  by  Banner  and  Van  Uitert^’^  for  crystal  pulling  by  the  Kyropoulos 
technique . 

The  mole  fractions  of  the  starting  oxides  were  PbO:0.6M,  MgO:0.2M,  Nb.,0^: 
0.06M,  B20y0.14M.  Care  must  be  taken  in  precalcining  the  Mg 0 powder  to  avoid 
hydration.  The  oxides  were  wet  mixed,  then  melted  down  into  a platinum  crucible. 

The  temperature  program  involved  a preheat  to  400°C,  held  at  that  tempera- 
ture for  several  hours  to  fully  oxidize  the  lead  and  reduce  the  likelihood  of 
reduced  PbO  reacting  with  the  crucible,  rapid  heating  to  1150°C  followed  by  a 
soak  of  2 hours  at  temperature.  The  melt  was  cooled  at  3°/hour  to  950°C,  at 
10°/hour  to  850°C,  followed  by  rapid  cooling  to  room  temperature. 

The  flux  was  not  poured  off  in  these  initial  runs,  and  the  crystals  were 
extracted  from  the  matrix  by  boiling  in  20%  nitric  acid  solution. 

The  crystals  grown  were  light  yellow  in  color  with  nearly  perfect  <100> 
cube  faces  and  edge  lengths  in  the  range  1-2  mm  (Fig.  5.1).  Transmitted  light 
microscopy  shows  most  of  the  crystals  to  be  free  from  major  flux  inclusions, 
isotropic  in  polarized  light  and  free  from  any  evidence  of  a gross  twin  structure. 
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In  the  thicker  samples,  the  edge  regions  exhibit  parallel  extinction, 
but  this  is  probably  due  to  total  internal  reflection  at  the  inside  surfaces. 
The  body  of  the  crystals  appears  strain  and  extinction  free. 

Preliminary  dielectric  measurements  show  a high  dielectric  maximum  at 
temperature  just  below  room  temperature  which  appears  dispersive  as  is  expected 
in  the  lead  magnesium  niobate. 

5 . 4 Growth  of  P2T  Solid  Solution  Crystal s 

It  should  be  stressed  at  the  outset  that  the  goal  of  this  longer  range 
program  is  not  to  produce  large  crystals  of  PZT  for  any  practical  application. 
The  requirement  is  for  very  small  homogeneous  crystals  with  simple  twin  struc- 
tures which  can  be  de-twinned  mechanically  and  electrically,  so  that  the  single 
domain  properties  are  accessible  for  study.  Thus  the  critical  features  are 
homogeneity  in  the  cation  distribution,  perfection  and  morphology.  Size  is  of 
very  minor  importance. 

Writing  the  formula  for  the  PZT  family  in  the  form  Pb2r  Ti,,  .0~,  the 

J x (1-x)  3 

values  of  x covering  the  range  0.10  to  0.80  are  of  most  importance,  since 
this  covers  the  full  range  of  simple  proper  ferroelectric  compositions  and 
encompasses  the  important  region  of  morphotropy. 

In  the  initial  phase  of  the  present  approach,  four  avenues  are  presently 
being  explored.  For  PZT  feed  stock,  both  simple  mixed  oxides  and  precalcined 
"homogenized"  PZTs  are  being  tried.  Two  flux  systems  are  also  being  employed. 
Tlie  first  comprises  a mixture  of  lead  fluoride  potassium  fluoride  and  lead 
phosphate,  which  has  been  used  successfully  by  Rosolowski  and  Arendt 
for  the  growth  of  x = 0.50  PZTs.  The  second  was  suggested  by  recent  work  at 
General  Electric  Company,  where  highly  reproducible  PZT  powders  are  being 
produced  from  a KCl:NaCl  eutectic  melt. 
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To  date,  only  some  10  growth  runs  have  been  attempted.  In  each  case 
a fixed  15:85  weight  ratio  of  PZT:flux  was  used  and  heating  was  carried  out 
in  a small  50  cc  platinum  crucible. 

For  the  fluoride:phosphate  flux  the  temperature  cycle  involved  rapid 
initial  heating  (3  to  3-1/2  hours),  a 3-hour  soak  at  1210°C,  cooling  to  810°C 
over  a period  of  1-1/2  to  5 hours,  then  a rapid  cool  down  to  room  temperature. 
In  several  runs  the  flux  was  poured  off  at  810°C;  for  one  or  two,  the  flux 
was  solidified  in  the  crucible. 

In  each  case,  small  single  crystals  were  obtained.  A typical  example 
is  shown  in  Fig.  5.2  from  a run  with  x = 0.80.  Here  the  optical  properties 
show  clear  straight  extinction  in  polarized  light,  indicating  that  the  grown 
composition  is  in  the  rhombohedral  phase  field  at  room  temperature.  High 
temperature  microscopy  is  now  underway  with  these  crystals  to  determine  the 
homogeneity  of  the  Curie  temperature  within  a single  crystal,  and  between 
crystals  of  the  same  batch. 

Further  effort  is  obviously  required  to  modify  the  growth  process  to 
improve  the  perfection. 

For  the  KCl:NaCl  flux,  less  work  has  been  carried  out,  and  the  initial 
results  do  not  look  encouraging.  Soak  tests  carried  out  on  the  mixed  oxide 
compositions  at  950°C  suggest  that  while  the  PbO  is  rapidly  digested  in  the 
flux,  ZrC^  and  TiO,  remain.  We  propose  now  to  experiment  with  organically 
prepared  PZTs  as  the  starting  material  to  see  if  this  will  dissolve  in  the 
flux.  If  these  experiments  are  not  successful,  all  effort  will  be  concentrated 
on  the  other  flux  system. 

5.5  Reference 


5.1  W.A.  Bonner,  L.G.  Van  Uitert.  Mat.  Res.  Bull.  2,  131  (1967). 


Phonbohedral  crystal  from  Run  Rol  -f-d2-10?  shown  in  a)  transmitted 
unpolarised  linht  and  h)  transmitted  polarised  liaht  to 
illustrate  the  development  of  domains. 
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6.  STUDIES  OF  FERROELECTRIC  BICRYSTALS 

6 . 1 Introduction 

The  initial  objective  of  these  studies  is  to  produce  "joining"  between 
differently  oriented  ferroelectric  crystals,  along  an  extended  planar  boundary, 
which  will  simulate  as  far  as  possible  the  character  of  a single  grain  boundary 
in  a ferroelectric  ceramic.  Such  planar  macro  boundaries  would  permit  the 
measurement  of  dielectric,  optical  and  elastic  properties  and  a possibility 
for  untangling  the  complex  interaction  of  elastic  and  electric  boundary  con- 
ditions on  the  single  crystal  properties,  and  on  the  evolution  of  the  domain 
structures  in  the  ferroelectric  phases. 

In  normal  ferroelectric  ceramics  the  grain  boundaries  do  encompass  quite 
a range  of  perfection,  of  width,  and  of  chemical  character  from  the  near  perfect 
co-joining  of  the  narrow  angle  grain  boundary  to  the  very  broad  second  phase 
inclusion  in  some  deliberately  diphasic  systems.  It  is  thus  of  interest  to 
look  at  a range  of  techniques  and  of  ferroelectric  systems  before  making  de- 
tailed choices  for  specific  study. 

6 . 2 Experimental  Approach 

Two  possible  techniques  for  producing  planar  grain  boundaries  have  already 

(6.1) 

been  proposed  viz:  direct  Czochralski  growth  of  two  oriented  seed  crystals 

in  a suitable  melt  until  the  crystals  impinge,  or  alternatively  diffusion 
bonding  between  two  oriented  single  crystals  under  pressure  at  high  temperature. 
The  two  methods  are  compared  in  Table  6.1.  In  view  of  the  comparative  simpli- 
city of  the  diffusion  bonding  technique,  this  approach  has  been  attempted  first. 

(6.2) 

It  has  been  reported  that  in  the  presence  of  a suitable  solvent  many 
powdered  materials  can  be  pressed  to  high  density  at  temperatures  at  or  below 
room  temperature.  Examples  include  the  ice/methanol  and  sodium  chloride/water 
systems.  The  solvent  acts  as  a high  dlffusivity  medium  for  material  transfer 
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Table  6.1  Comparison  of  Bicrystal  Preparation  Methods 


Method 


1.  Czochralski 


2.  Diffusion 
bonding 


Requires  accurately 
oriented,  flaw-free  seed 
crystals,  a suitable 
melt,  and  development 
of  the  correct  growth 
techniques  for  each  com- 
position required. 


Can  be  prepared  from  a 
single  crystal  by  cutting 
the  crystal  and  grinding 
to  give  the  required 
orientation.  Requires  a 
material  with  sufficient 
atomic  mobility  at  the 
bonding  temperature. 


Product 


Grain  boundary  free  from 
secondary  boundaries  intro- 
duced by  recrystallization. 
Possible  contamination  from 
the  melt  if  a flux  is  used. 

No  direct  comparison  between 
the  properties  of  the  orig- 
inal crystals  and  the  bi- 
crystal is  possible. 

Grain  boundary  is  directly 
comparable  to  those  produced 
in  sintering  and  hot  press- 
ing. Second  phases  can  be 
incorporated  in  the  boundary 
prior  to  diffusion  bonding. 
Direct  comparison  between 
the  properties  of  the  orig- 
inal crystal  and  the  bicrystal 
can  be  made.  Contamination 
during  grinding  and  polishing 
and  recrystallization  due  to 
strain  introduced  during  grind- 
ing may  occur. 
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during  deformation.  In  many  ways  the  process  is  analogous  to  liquid-phase 
sintering.  It  should  be  possible,  therefore,  to  prepare  bicrystals  from  water 
soluble  ferroics  such  as  TGS  and  KDP . Attempts  have  been  made  to  produce  TGS 
bicrystals  using  this  method. 

6.2.1  Diffusion  Bonding 

Specimens  for  diffusion  bonding  are  prepared  as  follows:  After  measuring 
its  dielectric  properties,  a single  crystal  is  cut  in  half  using  a string  saw 
to  minimize  structural  damage.  Thu  halves  are  then  ground  and  polished  on  one 
face  in  the  appropriate  orientations  until  surface  irregularities  are  reduced 
to  the  submicron  level.  Diffusion  distances  required  to  remove  them  entirely 
and  permit  bonding  are  then  similar  to  those  encountered  in  normal  sintering 
and  hot  pressing. 

A simple  hot  press,  shown  in  Fig.  6.1,  has  been  built  to  bond  the  crystals. 
The  specimen  is  heated  using  a platinum  40 % rhodium  wound  tube  furnace  capable 
of  attaining  1500°C.  Pressure  is  applied  by  deadweight  loading  through  two 
accurately  aligned  alumina  rods  with  parallel  ground  contact  faces.  Reaction 
between  the  alumina  rods  and  the  specimen  is  minimized  using  platinum  foil, 
which  also  serves  to  smooth  out  stress  gradients.  Materials  subject  to  com- 
position change  by  vaporization,  such  as  PZT,  are  sealed  in  platinum  foil  prior 

to  heating.  The  press  has  been  designed  to  apply  typical  hot-pressing  pressures 
-2  2 

(10-50  MNm  ) on  crystals  in  the  1-4  mm  size  range.  These  small  crystals  are 
cheap  and  comparatively  easily  grown  and  quite  flaw  free. 

The  hot  press  has  been  used  to  prepare  a bicrystal  of  quartz  at  1350°C 

_2 

and  10  MNm  pressure.  Bicrystals  of  BaTiO^,  SrTiO^,  KTN,  PZT  and  other  mate- 
rials will  be  made  as  crystals  become  available  from  flux  growth  work. 


6.2.2  Solution  Bonding 

Attempts  have  been  made  to  form  TGS  bicrystals  at  room  temperature  by 
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pressing  together  two  single  crystals  with  a layer  oi  saturated  TGS  solution 

between  them.  Tin  ml  ' th«  wo  :rystals  wei  ■ ! trst  highly  poli  hed 

using  powdered  tin  oxide  and  ethanol  as  a polishing  medium.  Pressures  were 
_o 

kept  below  1 MNm  “ since  higher  pressures  invariably  fractured  the  crystals. 

No  attempt  was  made  to  control  the  pressing  temperature.  After  allowing  l 
days  to  reach  equilibrium  the  crystals  were  desiccated  in  high  vacuum  (10 
at  60°C  to  evaporate  any  remaining  solution. 

So  far  this  "diffusion  bonding"  method  has  failed  to  produce  a strong  bond 
between  the  crystals.  On  exami  it  n,  evidence  of  severe  etching  was  found 
the  contact  faces  of  the  separated  crystals,  hu;  no  regions  where  bonding  had 
apparently  taken  place.  Accordingly,  a "crystal  growth"  approach  has  been 
adopted.  The  cryst  bed  previously  and  placed  togethei 

horizontally  with  the  thickest  possible  1 solution 

them,  but  without  pres  ter  two  days  at  room  tempi  ituri 

were  desiccated  under  va< u nn  overnight,  »ith  water  ai 

examined  by  scanning 

occurred  between  the  crystals  in  isolated  regions,  but  many  large  pores  remain. 
The  bonded  region  is  qud l 

probably  not  as  crystallographies:  v perfect  ar  the  original  crystals.  This 
may  be  associated  with  surface  damage  incurred  during  polish  ng. 

In  order  to  maximize  the  bonded  area  between  the  crystals  it  is  necessary 
to  raise  the  initial  concentration  of  tht  TGS  solution  as  much  as  possible. 
Unfortunately,  a high  solution  concent rat  i c precludes  etching  of  the  TGS 
crystals  and  prevents  the  formation  of  cK  :n  contact  surfaces.  Using  the 
simple  model  in  Fig.  6.3a,  iL  can  be  shown  oat  for  unit  area  of  TGS/solution 
interface , 
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FIGURE  6-3 (b)  (From  reference  6-A) 
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where  AL  is  the  thickness  of  TGS  crystal  removed  by  etching,  AC  is  the  change 
in  TGS  solution  concentration,  L is  the  initial  solution  thickness  and  p is  the 
density  of  TGS.  AL  can  be  increased  by  raising  L (which  can  be  accomplished  by 
making  the  upper  crystal  as  thin  as  possible  to  minimize  pressure  during  bonding), 
or  by  raising  AC. 

Figure  6.3b  shows  the  increase  in  TGS  solubility  with  solution  temperature. 
The  graph  indicates  that  raising  the  growth  temperature  of  the  bicrystal  would 
permit  the  maximum  room  temperature  solution  concentration,  and  yet  give  high  AC, 
and  hence  AL,  on  heating.  At  equilibrium  the  quantity  of  TGS  in  solution  between 
the  crystals  would  be  greatly  increased  compared  to  the  room  temperature  concen- 
tration, and  a larger  bonded  area  would  result  on  desiccation.  According  to 
Fig.  6.3b  the  growth  temperature  must  be  as  high  as  possible  to  obtain  maximum 
concentration  change.  However,  there  is  some  evidence11  ' that  TGS  decomposes 
slowly  above  60°C  and  therefore  this  temperature  probably  represents  the  practical 
upper  limit.  Experiments  to  test  some  of  these  ideas  are  at  present  in  progress. 

6.2.3  Characterization  of  Real  and  Artificial  Grain  Boundaries 

In  order  to  ascertain  the  types  of  boundaries  which  are  geing  generated 
so  as  to  be  able  to  compare  them  with  the  appropriate  class  of  ceramic 
grain  boundaries  in  polycrystalline  samples,  it  is  necessary  to  have  the 
capability  to  examine  the  boundary  region  on  an  atomic  scale.  Transmission 
electron  microscopy  is  the  only  technique  capable  of  providing  the  resolution 
required . 

To  facilitate  these  studies  a cooperative  program  has  been  initiated  with 
the  Department  of  Materials  Science  at  the  University  of  Virginia  in  Charlottes- 
ville. The  Department  in  Virginia  has  excellent  transmission  microscope  facili- 
ties under  the  direction  of  the  Department  Chairman,  Dr.  K.  Lawless,  and  is 
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concerned  with  developing  a regional  center  of  excellence  in  this  area. 

To  initiate  cooperation,  a study  of  the  grain  boundaries  in  PZT  ceramics 
was  begun  in  October  1978,  centered  on  the  question  of  domain  wall  continuity 
through  the  rather  broad  second-phase  boundary  material  which  is  evident  in 
SEM  studies  of  this  composition. 

We  believe  that  the  commonality  of  interests  between  the  faculty  in 
Virginia  who  wish  to  work  more  in  the  ceramic  systems,  and  have  the  tools  and 
expertise  for  excellent  transmission  work,  and  from  our  own  understanding  of 
the  electroceramic  problems  and  urgent  need  for  the  high-resolution  studies, 
an  excellent  joint  program  can  be  carried  forward. 
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7.  PHENOMENOLOGICAL  STUDIES 
7 . 1 Introduction 

In  simple  proper  ferroelectric  crystals,  Landau:Devonshire:Ginsburgh 
phenomenological  theory  provides  a simple  framework  which  allows  the  thermo- 
dynamic interrelations  between  elastic,  dielectric,  and  thermal  properties  of 
the  ferroelectric  species  accessible  from  any  paraelectric  prototype  structure 
to  be  related  to  the  first  and  higher  order  stiffnesses  and  coupling  parameters 
of  the  prototype.  From  the  Taylor  series  expansion  of  the  free  energy  the 
first  partial  derivatives  can  be  used  to  determine  possible  stable  ferroic 
species,  while  the  second  partial  derivatives  interrelate  the  properties  of 
species  and  prototype.  Unfortunately,  however,  the  sixth  order  polynomial  forms 
which  are  required  to  describe  the  weak  first  order  transitions  characteristic 
of  many  ferroic  species  lead  to  rather  cumbersome  polynomial  equations  which 
become  very  tedious  to  explore,  both  for  phase  stability  and  property  tensors. 

When  considering  the  fact  that  many  ferroelectrics  of  practical  interest 
are  selected  from  solid  solutions  of  end  member  components  which  are  proper 
ferroelectrics,  and  that  in  a solid  solution  family,  the  coefficients  in  the 
expansion  of  the  free  energy  mutate  only  slowly  and  continuously  with  composi- 
tion change,  it  is  evident  that  if  the  mathematical  manipulation  could  be 
automated,  the  phenomenological  approach  could  provide  an  effective  systematic 
manner  for  interpolating  in  the  exploratory  phase  diagrams  for  new  complex 
multicomponent  solid  solution  systems  and  for  predicting  the  single  domain  prop- 
erties of  the  ferroic  species  and  their  changes  with  changing  elastic,  dielectric 
and  thermal  boundary  conditions. 

Over  the  past  three  years,  we  have  expended  a considerable  effort  to 
explore  computational  techniques  whereby  the  tedium  of  the  method  could  be 
reduced  and  thus  its  potential  utility  enhanced.  Using  the  ADAGE  computer 
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graphics  system  at  PEnn  State,  programming  has  been  developed  for  the  Devon- 
shire form  of  the  elastic  Gibbs  function  for  the  perovskite  structure  ferro- 
electrics  whose  prototype  is  in  point  group  m3m.  The  ADAGE  which  is  a mixed 
analog  digital  system  has  been  set  up  so  that  the  dielectric  stiffness  and 
higher  order  stiffness  coefficients  and  their  temperature  coefficients  are 
under  dial  control.  For  any  choice  of  these  coefficients,  the  program  will 
map  the  energies  of  possible  stable  states  as  a function  of  temperature,  and 
calculate  polarizations  and  permittivities  in  the  resulting  possible  stable 
species. 

To  aid  in  the  exploration  of  solid  solutions,  coefficients  which  reproduce 
the  stable  phases  and  properties  of  the  end  members  can  be  set  up,  then  the 
machine  can  be  arranged  to  increment  the  coefficients  between  these  end  member 
values  plotting  both  AG:T  and  effective  phase  composition  vs  T diagrams. 

The  focus  of  this  work  has  been  to  develop  an  energy  function  which  would 
describe  the  single  domain  properties  of  all  simple  f erroelectrics  in  the  pure 
PbZrO^ : PbTiO^  solid  solution  system,  and  this  task  is  now  virtually  completed. 
The  chosen  function  reproduces  the  unusual  feature  of  morphotropy  between 
tetragonal  and  rhombohedral  phases,  reproduces  the  observed  spontaneous  strains 
from  realistic  values  of  the  electrostrict ion  constants  which  are  almost  inde- 
pendent of  composition,  gives  reasonable  predictions  of  the  single  domain 
dielectric  and  piezoelectric  constants  and  magnitudes  of  the  free  energies  in 
the  different  phases  which  agree  well  with  those  derived  from  thermal  measure- 
ments . 

7 . 2 Foci  of  Present  Phenomenological  Studies 
7.2.1  Perovskite  Ferroelectrics 

(a)  Using  the  coefficients  determined  in  the  earlier  study  for  composi- 


tions very  close  to  the  morphotropic  boundary,  the  effects  of  changing  the 
elastic  boundary  conditions  are  being  explored. 
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(b)  For  these  morphotropic  compositions  the  polarization  surfaces  of 
constant  elastic  Gibbs  free  energy  are  being  determined  as  a function  of  tem- 
perature and  of  composition.  The  objective  of  this  study  is  to  determine  the 
trajectory  of  the  polarization  which  gives  the  lowest  energy  switching  from 
rhombohedral  to  tetragonal  states. 

(c)  Modifications  of  the  Curie  temperature  and  higher  order  stiffnesses 
are  being  studied  which  mimic  the  behavior  of  the  PLZT  family. 

7.2.2  Tungsten  Bronze  and  Bismuth  Oxide  Layer  Structures 

For  the  very  large  families  of  tungsten  bronze  and  bismuth  oxide  layer 
structure  ferroelectrics,  the  prototype  symmetry  is  4/mmm.  In  both  families, 
it  appears  that  there  are  simple  proper  ferroelectrics  and  also  improper  ferroic 
species.  Thus  the  elastic  Gibbs  function  is  necessarily  more  complex.  Formally 
it  may  be  expressed  in  the  following  manner: 

AG  = 3^  (np)  + <I>2(XP)  + 3>3(XPrip) 

where  3>^  is  a function  of  the  independent  components  HP  of  the  order  parameters 
which  characterize  the  improper  ferroic  phase. 

$2  contains  the  contributions  related  to  the  elastic,  dielectric,  and 
elasto-electric  (electrostrictive)  energies  of  the  crystal. 

represents  the  coupling  between  the  order  parameters  for  the  improper 
coordinates  and  the  stress  and  polarization  components. 

The  form  chosen  for  the  function  3^  will  depend  on  the  nature  of  the  im- 
proper phases  displayed  by  the  system.  In  barium  sodium  niobate,  which  exhibits 
an  improper  ferroelastic  transition  at  543  K,  it  would  have  the  form 

= 1/2  a(n2+P2)  + 1/4  e1(n4+P4)  + 1/2  e2n2P2 

+ 1/2  Y^r^+P6)  + 1/2  Y2n2p2(n2+p2) 
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The  function  $2  expanded  to  include  the  first  sixth  order  terms  has  the 
general  form: 

*2  " “l  <Pf+P2>  + “3  P3  + ail  (PJ+P2>  + a33P3 


+ aX  (P^p^  + p^p^\  + aX  + nX 

ai3  (P1P3  + P2P?  + a!2PlP2  + a333P3 


+ “ni<Pl+P2)  - 1/2  S11<X1+X2>  - S12X1X2 


- SP3(X1+X2)X3  - 1/2  SP3X23  - 1/4  SP4(X2-fX25) 


+ Q13(pfx3+p*x3)  + Q31(P3X1+p2x2) 


+ Q33P3X3  + WP2P3X4+  V3V 


+ Q66  P1P2X6  ' 


The  function  4>3  would  again  depend  on  the  nature  of  the  impropriety,  and  for 
barium  sodium  niobate  the  lowest  order  interaction  terms  would  give 

*3  = ei  ^X1+x2) (n2+p2)  + e2x3(n2+P2) 


2222  222 

+ Pj_  (p3+p2)(h  +p  ) + p2  p3  (n  +p  ) 


+ v3  (x3-x2) (n2-P2)  + v2x6np 


At  any  temperature  the  stable  states  of  the  crystal  are  given  from  the 


conditions 
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SAG  = SAC 
Sri  Sp 


S2AG  _ S2AG 

2 2 > U 
Sn  Sp 


s2ag  s2ag  s2ag 

Sn2  Sp2  3n3p~ 


SAG  _ SAG 

3X.  SP . 

1 j 

X X 

It  is  expected  that  the  coefficients  cx^,  and  a may  be  strongly  and 
linearly  temperature  dependent  with  the  forms 


al  “ aUo)  (T"01) 


a3  a3(o)  (T  °2) 


a = a (T-0.) 
o 3 


Depending  on  the  magnitudes  of  0^,  02  and  0^  different  ferroelectric  and 
improper  ferroelastic  species  are  accessible  for  the  system. 

The  present  task  on  this  program  is  concerned  primarily  with  the  function 
$2  and  is  limited  to  the  description  of  simple  proper  ferroelectric  phases. 

The  material  family  of  most  interest  is  Pbfrt^O^,  Pb2NaNb^0^^  and  the  solid 
solutions  between  Pb2NaNb^0^^  and  Ba2NaNb^0^,_ . 


Data  are  being  assembled  and  the  programming  developed  to  describe  this 
solid  solution  system,  using  the  same  methods  applied  earlier  to  the  PZT  family. 
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8.  PROGRAM  ORGANIZATION:  PERSONNEL  AND  EQUIPMENT 
8.1  Personnel 

The  magnitude  and  scope  of  the  present  effort  have  required  a rapid 
build-up  of  personnel  and  equipment  for  the  new  topics  of  study  which  are 
being  undertaken.  Over  this  initial  period  we  have  been  most  fortunate  to  be 
able  to  draw  upon  the  services  of  two  distinguished  visitors,  K.  Uchino  from 
the  Tokyo  Institute  of  Technology,  Japan,  and  A.  Nicol  from  Birmingham  Univer- 
sity in  England.  Uchino  has  been  responsible  for  the  rapid  assembly  and  testing 
of  the  new  AC  precision  dilatometer.  Nicol,  who  spent  three  months  in  MRL  over 
the  summer,  advised  upon  and  devised  the  new  temperature  control  equipment  for 
crystal  growth,  and  helped  initiate  the  PZT  growth  studies. 

For  the  theoretical  work  on  electrostri ction.  Professor  Barsch  has  been 
joined  by  Dr.  N.  Achar.  Achar  has  had  considerable  earlier  experience  in 
problems  involving  fundamental  lattice  calculations,  and  has  worked  closely 
with  liarsch  on  several  earlier  theoretical  studies. 

Work  on  the  more  practical  electrostrictors  is  being  carried  forward  by 
two  graduate  assistants:  N.  Setter,  who  has  been  growing  Pb^lft^MgOg  single 
crystals  and  studying  order: disorder  phenomena  in  complex  perovskites,  and 
S.J.  Jang,  who  is  completing  studies  of  electrostriction  in  the  Pb^MgNb^O^: 
PbTiO^  ceramic  system. 

Overall  responsibility  for  all  electrostriction  work  is  with  Professor 
Cross,  who  is  working  closely  with  Uchino  and  Barsch  to  coordinate  experimental 
and  theoretical  approaches. 

The  diphasic  studies  under  the  direction  of  Professor  Newnham  have  been 
joined  by  Dr.  B.  Hardiman,  a new  faculty  member  who  has  joined  MRL  from  the 
A.G.  Telefunken  Laboratories  in  Hamburg.  Dr.  Hardiman  is  working  with  newer 
replication  techniques  to  produce  composites  of  controlled  connectivity  from 
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PZT  composites  are  being  carried  forward  by  Mr.  K.  Klicker,  who  has  just 
completed  his  M.S.  in  Solid  State  Science. 

Newnham's  work  on  grain  oriented  layer  structure  ceramics  is  also  assisted 
by  Mr.  M.  Holmes,  a graduate  assistant  working  on  a Bismuth  Institute  fellow- 
ship at  MRL. 

Crystal  growth  studies,  which  were  initiated  by  Dr.  Nicol  are  now  being 
carried  forward  by  Ms.  Susan  Stephenson,  a full-time  research  aide.  Susan  will 
be  the  back-up  person  for  all  growth  work  and  is  cooperating  with  Dr.  M.  Brun, 
who  is  the  professional  crystal  grower  presently  in  residence  at  MRL. 

Bicrystal  developments  are  being  pushed  by  Dr.  L.  Bowen,  a postdoctoral 
fellow  from  Leeds  University.  Dr.  Bowen  is  being  assisted  by  Susan  Stephenson 
in  the  water  soluble  regrowth  studies.  This  work  is  supervised  by  Dr.  J.V. 
Biggers. 

Phenomenological  work  on  this  program  is  being  done  by  Mr.  A.  Amin,  a 
graduate  assistant.  He  is  cooperating  closely  with  Dr.  H.  McKinstry  and  Dr. 

R.  Betsch,  who  are  developing  phenomenological  methods  for  understanding  the 
elastic  and  thermal  expansion  properties  of  tungsten  bronze  type  f erroelectrics 
on  other  program  funds. 

Dielectric  and  piezoelectric  measurements  are  coordinated  for  the  whole 
program  by  Dr.  W.  Schulze,  and  all  ceramic  preparative  and  characterization 
studies  are  under  the  direction  of  Dr.  J.V.  Biggers. 

8 . 2 Equipment 

The  major  investment  in  new  equipment  has  been  concerned  with  modernization 
and  automation  of  the  dielectric  and  piezoelectric  measurements.  After  a 
careful  review  of  competing  DEC  and  hobby-type  digital  equipment,  it  was 
decided  to  base  the  main  control  around  an  HP  model  9825  system.  Unfortunately, 
the  delivery  schedules  for  the  new  system  and  interface  equipment  were  such 
that  it  was  not  possible  to  put  them  into  operation  over  the  first  year  of 
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the  contract.  All  major  items  including  the  calculator,  plotter,  printer, 
programmer,  and  associated  accessories  are  now  scheduled  for  delivery  and 
we  anticipate  automating  our  dielectric  measurement  capability  early  in  the 
new  year. 

Three  other  major  additions  have  been  made  over  the  previous  year. 

1.  Major  equipment  required  for  the  new  AC  precision  dilatometer  described 
in  Section  2 have  been  installed. 

2.  The  complete  y data  track  program  control  arrangement  from  Research 
Industries  has  been  installed  and  is  operational. 

3.  In  cooperation  with  the  MRL  cement  group  under  Dr.  D.  Roy,  a Haake 
model  RV-3  rotovisco-viscosimeter  with  a wide  range  of  accessories  for  vis- 
cosities in  the  range  of  interest  has  been  purchased  and  installed. 

With  the  very  comprehensive  preparation  and  characterization  equipment 
now  available  in  the  laboratory,  the  group  is  in  excellent  shape  for  the 
developing  program.  For  the  more  detailed  study  of  resonant  devices  we  anti- 
cipate the  need  for  an  improved  network  analyzer,  and  this  will  be  a major 
item  required  in  the  new  year. 

8 . 3 Advisory  Committee 

The  first  meeting  of  the  ONR  Advisory  Committee  for  this  program  was  held 
at  Penn  State  on  Wednesday,  27th  September  1978.  Members  present  were 
A.  Glass,  Bell  Telephone  Laboratories;  H.  Anderson,  University  of  Missouri- 
Rolla;  B.  Jaffe,  Vernitron  Corporation,  and  V.  Smith,  Philips  Laboratories 
(for  S.K.  Kurtz).  Unable  to  attend  were  G.  Kino,  Stanford  University  and 
D.  Berlincourt,  Channel  Products.  Navy  representatives  at  the  meeting  were 
R.  Pohanka,  ONR  and  NRL,  Washington;  P.  Smith,  NRL,  Washington;  C.  LeBlanc, 

NUSC,  New  London,  and  R.  Wollett,  NUSC,  New  London.  Guests  present  at  the 
meeting  were  K.  Okazaki,  National  Defense  Academy,  Tokyo;  G.  Taylor,  Rutgers 
University  (Editor,  Ferroelectrics) , and  D.  Buckner,  Gulton  Industries. 
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The  Committee  were  given  an  overview  of  the  new  program  during  the 
morning  of  the  27th  of  September.  The  afternoon  was  spent  in  discussions 
and  in  exchange  of  information  concerning  the  development  of  f erroelectrics 
and  piezoelectrics  in  the  USA  and  Japan.  These  proceedings  have  been  reported 


separately. 


